Field-Induced Breakup of Emulsion Droplets Stabilized by Colloidal Particles 
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We simulate the response of a particle- stabilized emulsion droplet in an external force field, such 
as gravity, acting equally on all N particles. We show that the field strength required for breakup (at 
fixed initial area fraction) decreases markedly with droplet size, because the forces act cumulatively, 
not individually, to detach the interfacial particles. The breakup mode involves the collective desta- 
bilization of a solidified particle raft occupying the lower part of the droplet, leading to a critical 
force per particle that scales approximately as iV -1 ' 2 . 

PACS numbers: 82.70.Kj (emulsions); 82.70.Dd (colloids); 81.40.Np (failure of materials) 
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The stabilization of emulsions is an important issue in 
industry, where it controls the functionality and shelf-life 
of many products [l|. Alongside traditional surfactants, 
colloidal particles or nanoparticles have long been used to 
stabilize emulsion droplets (Pickering emulsions) [2[ , with 
several emerging applications that exploit the nonequilib- 
rium nature of their adsorption to the fluid- fluid interface 
[3|45|. However this also poses obstacles to formulation: 
to maintain emulsion stability, it is essential that the col- 
loidal particles remain sequestered on the interface, even 
when the system is subject to external challenges [3, 6]. 

Brownian motion is generally insufficient to cause de- 
tachment [SlZI, so we ignore it here. In contrast, external 
body forces acting equally on each particle can, if strong 
enough, overcome the capillary forces holding them in 
place. Destabilizing forces can be caused by magnetic or 
electric fields or their gradients, and also by gravity. For 
simplicity we refer mainly to gravity in what follows. 

A groundbreaking experimental study was performed 
by Melle et al. [8[, using decane droplets in water sta- 
bilized by paramagnetic particles. They demonstrated, 
using a simple bar magnet, controlled field-gradient in- 
duced breakup of an otherwise highly stable emulsion. 
A second study uses instead the dielectrophoretic force 
caused by an electric field gradient [9|. A third recent 
study has found complex dynamics in droplets stabilized 
by supra-colloidal (400/xm) particles under gravity [10j . 
(Scaling by the Bond number, defined below, indicates 
that the same physics should arise for, e.g., 1/im particles 
under bench-top centrifugation.) Although field- induced 
breakup is simultaneously a threat to an emulsion's sta- 
bility and a promising route to its controlled breakage, 
the dependence of the critical field for breakup on the 
droplet size is far from understood. 

Previously we performed lattice Boltzmann (LB) sim- 
ulations on emulsion structures (primarily bicontinuous 
|ZL LUL| ) stabilized by magnetic colloids in the presence of 
a magnetic field gradient [12, LL3| . There we argued that 



the effect of the force does not act separately on each 
particle but is cumulative across the structure. (This ex- 
plains the very low field-gradients required in [8].) By 
that argument, a major fraction of the weight of all par- 
ticles could be loaded onto a single 'keystone' particle at 
the droplet bottom which should then detach, initiating 
breakup. In this paper, we study in detail the breakup 
dynamics of isolated emulsion droplets coated with N (es- 
sentially hard-sphere) particles under gravity. We find 
a scaling of the critical force with N that is inconsis- 
tent with the 'keystone' argument, and associated instead 
with a failure mode in which a raft of jammed particles 
detaches from the droplet base taking liquid with it. 

Control parameters: For a single particle adsorbed 
onto a flat horizontal interface between two fluids of equal 
density, detachment is controlled by the Bond number 
[14j Bo = 3F/(47rcra) where a is the particle radius, a 
the fluid-fluid interfacial tension, and F the normal de- 
tachment force (in gravity, F is the buoyancy force on 
the particle). Detachment of an isolated particle (N = 1) 
occurs at a critical Bond number Boc(l) which depends 
on the fluid-solid interfacial tensions through the contact 
angle 6. In this work we consider only neutral wetting 
(0 = 90°) for which Bo c (l) = 3/4 pjHHj. The critical 
Bond number for detachment of one particle from a fluid 
sphere of radius R should also depend on r = R/a; to 
allow for this we define Boc(l, r) such that Boc(l, oo) = 
Bo c (l). The numerically determined Bo c (l,r) addition- 
ally reflects various sources of discretization error (e.g., 
the fact that in LB the particle radius is not much larger 
than the interfacial width |17j). 

Simulations: The LB methodology used in this work 



is by now standard [12|, LL5|, LL8| and summarized in [1 



The initial condition in our simulations comprises a fluid 
droplet of radius i?, immersed in another fluid of equal 
mass density, with N neutrally wetting particles present 
on its surface. We choose a near-constant surface area 
fraction a = N7ra 2 /47rR 2 = N/4r 2 = 0.53 - 0.54; this 



is high enough to represent a stable Pickering emulsion 
droplet but low enough to allow particles to rearrange on 
application of even a weak body force. The chosen N val- 
ues are TV = 42, 50, 102, 162, 194 and 700. Each droplet 
is placed in a box, with periodic boundary conditions 
(a repeat with closed-box boundaries for N = 700 gave 
near-identical results), whose LB lattice size (depending 
on R) lies between 64 2 x 128 and 200 2 x 400. 

We have examined three different initial arrangements 
of particles: one disordered, one a regular hexagonal 
packing and one a regular square packing (subject to the 
defects required to tile these onto a sphere). For the 
ordered configurations, initial particle coordinates lie on 
the vertices of a suitable polyhedron. The particles and 
droplet are then relaxed in the absence of body forces 
(F = 0) to allow the interfacial structure to equilibrate. 
Once this is done, a vertical body force F is switched 
on; this acts directly on each particle. To prevent overall 
translation of the droplet's center of mass, a balancing 
buoyancy force NF is applied as a uniform upward force 
density on the fluid nodes occupied by the droplet inte- 
rior. The structure of the droplet is then evolved by the 
LB algorithm which faithfully includes the effects of fluid 
flow and hydrodynamic interactions. 

Results: We define the critical Bond number Bo c (iV, r) 
as the threshold beyond which one or more particles 
break away from a droplet covered by N particles. For Bo 
< Boc(N, r) the droplet instead achieves a final steady 
state that is distorted but not ruptured. Since we work at 
constant a = N/Ar 2 ~ 0.53, we now drop the argument 
r and refer simply to Boc(iV). In practice we find par- 
ticle ejection always precedes disconnection of a droplet 
into two pieces and therefore used the time of first de- 
tachment, td, as a measure of when breakup occurs. (As 
time unit we adopt the Stokes time ts = 67rr]a 2 /F.) 

Starting from an equilibrated spherical droplet and 
switching on the external forces, particles start to move 
downwards (while the buoyancy force holds the droplet 
up). Upper sections of the droplet become bare [8|, LLOj 
and a dense interfacial 'sediment' is created towards the 
bottom of the droplet. As the sediment builds, forces are 
transmitted to the particles beneath. The lower part of 
the droplet can then become highly distorted. In con- 
trast to droplets initiated from a disordered configura- 
tion, which maintain a prolate spheroidal shape through- 
out the stable regime, the ordered arrangements give al- 
most symmetric force patterns in which the most unsta- 
ble particle lies at the tip of a cusp that develops at the 
droplet bottom. For Bo < Boc(iV), such a cusped con- 
figuration can remain stable indefinitely. Fig. [1] shows 
snapshots for both stable and unstable cases at the time 
td- Shown in red are those particles whose interparticle 
separation h < 0.04a. This criterion identifies the main 
load-bearing particles in the system [17j . 

Although the body force changes only by two percent 
between the highest stable and lowest unstable value ex- 
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FIG. 1. (Color online.) (a-h): Snapshots at N = 194. (a) 
The ordered initial configuration; (b)-(d) are snapshots at 
t = 37t s evolved from (a) with Bo = 0.050, 0.051 and 0.052 
respectively. Stable states are shown in (b) Bo = 0.050 with 
N c = 80, and (c) Bo = 0.051 with N c = 78. (d) Unstable 
state found at Bo = 0.052 with N c = 121 at t = t d = 37t s . 
(e) The disordered initial configuration; (f)-(h) are snapshots 
at t = 30t s evolved from (e) with Bo = 0.050, 0.053 and 0.055 
respectively. Stable states are shown in (f) Bo = 0.050 with 
iV c = 74, and (g) Bo = 0.053 with N c = 80. (h) Unstable 
state found at Bo = 0.055 with N c = 122 at t = t d = 30t s . 
(i)-(l): Snapshots for N=700. (i) The disordered initial con- 
figuration; (j)-(l) are snapshots at t = 14.5£s from the con- 
figuration (i) with Bo = 0.030, 0.033 and 0.034 respectively. 
Stable states are shown in (j) Bo = 0.030 with N c = 193, and 
(k)Bo = 0.033 with N c = 315. (1) Unstable state found at 
Bo = 0.034 with N c = 342 at t = t d = 14.5t s . Particles 
separated by h < 0.04a from a neighbor are shown in red. 



amined in Fig.[TJ there is a clear difference in their config- 
urations at the detachment time td- In both cases, close- 
packed clusters of N c load-bearing particles are observed 
at the bottom of droplets; while the geometry of these 
regions depends on initial conditions as detailed above, 
in the last stable droplet configuration A^ c c± 0AN for the 
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TABLE I. The time td of first particle detachment and the 
size N c of the load-bearing cluster immediately prior to td for 
unstable droplets of N — 194 with ordered and disordered 
initial states. Superscript * denotes Boc(iV) (numerically, the 
smallest Bo found for which breakup occurs). 



ordered and disordered cases. In the unstable state, im- 
mediately prior to particle ejection, 7V C ~ 0.67V is rather 
larger but again nearly independent of the initial con- 
dition. In Table U we report N c (td) for several values 
of Bo just above the threshold, in the case of ordered 
(square) and disordered packings with N = 194. We 
find N c > N/2 in all unstable cases, with rather little 
dependence on Bo (once the critical value is exceeded). 
As expected, however, a stronger body force (larger Bo) 
causes the detachment time td to decrease markedly. 

Fig. [2] depicts the time evolution of 7V C in the cases with 
TV = 194. This shows the buildup of the sedimented 
cluster after switching on gravity, followed by a period 
of saturation. Beyond t = 13ts, stable droplets reach 
steady state with roughly constant 7V C , although disor- 
dered cases show fluctuations due to reorganization of 
the randomly packed 2D sediment. In the unstable cases, 
such reorganization is always seen during the plateau re- 
gion; following a series of these plastic events N c starts 
increasing steadily until breakup occurs. 
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FIG. 2. (Color online.) The evolving cluster size, N c (t) for the 
simulations shown in Fig.[T](a-h). Left: N c (t) for the ordered 
initial configuration: black circles, Bo = 0.050, blue triangles 
Bo = 0.051, red squares Bo = 0.052. Right: disordered initial 
configuration: black circles Bo = 0.050, blue triangles Bo = 
0.053, red squares Bo = 0.055. For unstable cases, the plots 
end upon particle ejection. 
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FIG. 3. (Color online.) Log-log plot of critical Bond number 
for a droplet covered with N particles, Boc(iV), vs N. Blue 
diamonds: square ordered; green squares: disordered; black 
circles: hexagonal ordered. Blue dot-dashed line: best fit of 
all data to Bo c (iV) = AN~ 1/2 (A = 0.71). Red dotted line: 
best fit to Bo c (A0 = AN B (A = 0.38; B = -0.37). Black 
dashed line: bound found from A = 3/2 (see text). 



Scaling of critical Bond number: The above mecha- 
nism entails asssembly of a close-packed loadbearing clus- 
ter, involving roughly half the particles, which either 
maintains itself indefinitely (Bo < Boq(N)) or, after a 
series of plastic rearrangements, finally gives way, caus- 
ing particle ejection and droplet rupture (Bo > Boc(iV)). 
The droplet prior to this event is distorted significantly; 
its shape, while dependent in detail on the initial condi- 
tion, is roughly independent of TV, and involves order-one 
departures from a sphere. There is no sign of extreme 
cylindrical elongation as can arise in shear flow 2QL J21|. 

Although breakup is preceded by particle ejection, 



with all parameters studied here we do not see a stream 
of particles being ejected from the lower tip of the droplet 
as might be expected were the sediment of loadbearing 
particles to behave as a fluid. Such an outcome would al- 
low the weight of all N particles to be borne by a single 
'keystone' particle at the bottom of the droplet. Ejec- 
tion of that particle, possibly nucleating an instability of 
the entire structure, could then be expected at a criti- 
cal Bond number scaling as Boc~Boc(l, r)/N [12[. For 
our nearly hard-sphere particles [17], this is not observed 
here [24]; we instead find a scaling closer to TV -1 / 2 (see 
Fig. [3]). This is because the cluster has formed a jammed 
structure in which forces are transmitted laterally as well 
as vertically. This semi-rigid assembly then becomes col- 
lectively unstable (see Movies at [17j). 

Given the destabilization mode observed, we now esti- 
mate Boc by addressing the force balance on the close- 
packed sediment, viewed as a solid cluster. (See J22|, |23| 
for related arguments.) This has weight FN C ~ FN/2 
and to a first approximation occupies the lower half of 
a droplet that is only moderately distorted. Accordingly 
at Bo = Bo c this weight is just balanced by the surface 
tension force acting vertically across the equator of the 
deformed droplet, 2naR where R ~ R/2 is the radius at 
the equator. Bearing in mind that the initial surface cov- 
erage (here held constant) obeys a = Na 2 /AR 2 = N/Ar 2 
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suggests a critical Bond number scaling as TV -1 / 2 which 
gives a satisfactory first account of the simulation data. 
Acknowledgments: EGK thanks Siegfried Dietrich for 
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TABLE II. Critical Bond numbers (Boc(iV) for a droplet 
coated with N particles at initial coverage a ~ 0.53, and 
Boc(l, r) for a single particle on a droplet of the same radius); 
size of the load-bearing cluster prior to breakup; detachment 
time; breakup parameter A; and shape parameter R/R. Ini- 
tial configurations S, square, H, hexagonal, D, disordered. 



we then have the estimate Boc(iV, r) ~ 3R/Na or, in- 
troducing a negotiable geometrical constant A ~ 0.75, 
Bo c (N, r) = X/VaN. The predicted scaling Bo c ~ 7V" 1/2 
is tested in Fig. [3] which plots all data for Boc against N 
on a log- log plot. The best fit (treating data for all initial 
conditions equally) has A = 0.52. 

Table [II] shows all the Bo c values found in this study 
alongside run- by-run values of A. The residual increasing 
trend of A with N is consistent with the fact that the pre- 
breakup droplet shape becomes more spherical for larger 
N. Indeed the quantity XR/R, with ttR 2 defined as the 
cross sectional area of the deformed droplet at the height 
of the uppermost load-bearing particle, has no residual 
trend (see Table HI]) [251 ]. A free fit to the exponent gives 
Bo c ~ N-°- 37 , or R/R ~ TV 013 . We have no explana- 
tion for this weak residual trend, but note that it cannot 
persist to very large N unless droplets become oblate. 
Excluding that outcome, a bound Bo c > 3/(2y/aN) at 
large TV is found by balancing the tensile force 2iraR in 
an undeformed spherical configuration against the weight 
FN/2 of its lower hemisphere. (Any prolate deformation 
reduces the tensile force, decreasing Bo c .) This bound is 
plotted in Fig.[3j extrapolating our TV -037 fit, saturation 
is expected only at N ~ 10 4 . 

Conclusion: We have investigated by LB simulations 
the destabilization of an emulsion droplet coated with 
hard spheres under a body force that acts equally on all 
the particles [8l4l0| . We confirmed that the force required 
for destabilization is far less than would be expected if 
the force were to cause independent detachment of the 
particles from the fluid- fluid interface [12|. However, it 
is also far greater than would be expected if the weight 
of all N particles were to accumulate as a net detach- 
ment force on a single 'keystone' particle at the bottom 
of the droplet. That picture, which implicitly considers 
the stabilizing particle layer to comprise a 2D fluid, must 
now be replaced by one in which the lower part of that 
layer has solidified. A simple geometrical analysis then 
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FIELD-INDUCED BREAKUP OF EMULSION 

DROPLETS STABILIZED BY COLLOIDAL 

PARTICLES: ADDITIONAL ONLINE MATERIAL 

SIMULATION METHOD 

All simulations are performed using the LB method 
on a D3Q19 lattice for a binary fluid, with the bounce- 
back-on links method used to handle momentum ex- 
change with solid particles [18], as extended to the 
binary fluid case in [15j. The fluid- fluid interfacial 
tension and width are determined by the parameters 
of Ginzburg-Landau free energy functional F[ip] = 
J dV (tt^ 2 + "f ^ 4 + § (V^) 2 ), where ip is a composi- 
tional order parameter. Following [12|, ll9j we express 
these parameters in lattice units (LU) as — A = B = 
0.0011832 and k = 0.0023664; the interfacial thick- 
ness and the nominal interfacial tension are deduced as 
f = 2.0 LU and a = 1.5776 x 10~ 3 LU respectively. The 
fluid viscosity is set at r\ = 0.1 LU. Lattice units are as de- 
fined in [12j but for our purposes are arbitrary: all that 
matters is the Bond number. Note that the Reynolds 
number, which controls the relative importance of inertia 
compared to viscosity, remains small throughout our sim- 
ulations which represent the Stokes flow limit [12]. The 
sedimentation velocity of one particle v s can be used to 
define a capillary number Ca = r]v s /a = 2Bo/9 which 
remains small in all our simulations. 

From LB benchmark tests with this parameter set, 
the numerical value of surface tension is determined as 
a = 1.51 x 10 -3 LU which is 4% less than the nomi- 
nal value (Tq. (Discrepancies of this scale are common 
in LB (l2|.) Monodisperse neutral- wetting particles are 
introduced with hydrodynamic radius ah = 3.9 LU. This 
choice is imposed by the numerical resources available; 
note however that ah is only twice the interfacial width £. 
For that reason, one expects some discrepancy between 
the values of Boc(l, r) seen in simulations and those aris- 
ing for a structureless interface (£ — »> 0). The values 
tabulated in the main text are those found by numeri- 
cal simulations on a single particle residing on a droplet 
of radius R = ra. We have however confirmed that for 
a flat interface (r — >> oo), our methods recover the ex- 
pected result: indeed we find Boc(l) = 0.74 ~ 3/4. Par- 
ticle ejection is defined by the surface- average value of ip 
becoming sufficiently different from zero. 

Interparticle forces are mediated by a soft-core poten- 
tial U sc (h) chosen to prevent the surface-to-surface in- 
terparticle separation ft from becoming too small. (To 
model capillary forces correctly, it is important that 
the lattice Boltzmann nodes surrounding each parti- 
cle mainly represent fluid not solid sites.) We choose 
U sc (h) = u{h) - u{h c ) - (ft - h c )(j%) h=hc , where u{h) = 
7(fto/ft), and ft c is the maximum range of the interac- 
tion; we select 7 = 10 -4 , fto = 0.1, and h c = 0.25 [12j . 
In practice, the particles behave as nearly impenetrable 



spheres of radius a ~ ah — 4 LU. Finally, to allow for the 
discretization error in the hydrodynamic forces at small 
separations, a lubrication force (proportional to the nor- 
mal component of the interparticle velocity) acts when- 
ever h < 0.3 LU. The criterion ft < 0.15 LU = 0.04a is 
used to identify load-bearing particle contacts; this lies 
within the range of the soft-core potential defined above. 



ADDITIONAL RESULTS 

Figure H] shows snapshots of the particle configurations 
in much smaller droplets (N = 50), prepared initially 
in ordered (square lattice) and disordered states, at two 
different later times. These are for the lowest unstable 
Bo (~ Boq(N)) so that both droplets will subsequently 
break. The left pictures in Figures HJ^a) and [4fb) are 
configurations where N c has reached the near-saturation 
plateau value (compare Fig. [2]) showing semi-crystalline 
clusters at the droplet bottom. The right pictures in 
Figures HJa) and[H(b) are the configurations immediately 
before the first particle detachment. In each case there 
is a neck close to the droplet equator, with load-bearing 
particles below the neck. Although N here is quite small, 
these configurations are remarkably similar to those seen 
in Fig. m for N = 194 and even TV = 700. 

Figure [5] shows the breakup mode for an emulsion 
droplet (TV = 194) in which the hard-core interaction 
radius of the particles is unchanged but the radius they 
occupy on the interface is halved. For a given gravita- 
tional force per particle, this does not alter the force bal- 
ance argument for the collective detachment of the lower 
hemisphere. However it does reduce four- fold the maxi- 
mum total force that can be transmitted to the bottom- 
most particle without this detaching individually from 
the interface. In this case, a 'keystone' type failure mode 
is seen, in which particles stream off the bottom of the 
droplet without entrainment of the enclosed fluid. Note 
that similar physics could arise for partially wetting par- 
ticles with a contact angle far from 90 degrees: in this 
case also, the radius of the interfacial disc covered by 
each particle is significantly less than its hard sphere in- 
teraction radius. This again reduces the collective force 
sustainable without individual detachment, and therefore 
lowers stability with respect to the keystone detachment 
mode without having the same effect on the collective 
mode (in which a naive force balance again equates the 
total mass in the lower hemisphere to 27rRcr, independent 
of contact angle). 



MOVIES 

The following online movies are available showing the 
time evolution of stable and unstable droplets for various 
sizes and initial conditions. 



Movie 1. N = 194, ordered initial condition, stable 
(parameters as in Fig. lb). 

Movie 2. N = 194, ordered initial condition, unstable 
(parameters as in Fig. Id). 

Movie 3. N = 194, disordered initial condition, stable 
(parameters as in Fig. If). 

Movie 4. N = 194, disordered initial condition, unsta- 
ble (parameters as in Fig.lh). 

Movie 5. N = 700, disordered initial condition, stable 
(parameters as in Fig.lj). 

Movie 6. N = 700, disordered initial condition, unsta- 
ble (parameters as in Fig. 11). 

Movie 7. N = 194, ordered initial condition, unstable 
to 'keystone' type failure mode (parameters in Fig. 5c). 




(a) 




(b) 



FIG. 4. Comparison of particle-packing between ordered and 
disordered initial configurations at different time steps, (a) 
Snapshots of a droplet at t = 89ts on the left and at ta = 
lOOts on the right from initial ordered configuration, (b) 
Snapshots of a droplet at t = 5 Its (left) and at td = 59£s 
(right) from disordered configuration. 
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(b) 





(d) 



FIG. 5. Snapshots for droplets at N = 194 with hard core in- 
teraction radius twice the capillary radius, showing keystone 
detachment mode in place of the collective breakup. (The 
capillary radius is used in the definition of the Bond number; 
the hard core radius is used for the visualization.) (a) The or- 
dered initial configuration; (b)-(d) are snapshots at t — 11 At s 
evolved from (a) with Bo = 0.048, 0.064 and 0.11 respectively. 
(b) Bo = 0.048 with N c = 157 and N f = 2, (c) Bo = 0.064 
with N c = 160 and N f = 8. (d) Bo = 0.11 with N c = 157 
and N f = 21. 



